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Abstract 

Using the seesaw mechanism, we construct a model for the light-neutrino Ma- 
jorana mass matrix which yields trimaximal lepton mixing together with maximal 
CP violation and maximal atmospheric- neutrino mixing. We demonstrate that, 
in our model, the light-neutrino mass matrix retains its form under the one-loop 
renormalization-group evolution. With our neutrino mass matrix, the absolute neu- 
trino mass scale is a function of \U e s\ and of the atmospheric mass-squared difference. 
We study the effective mass in neutrinoless (3(3 decay as a function of | U e 3 | , showing 
that it contains a fourfold ambiguity. 
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1 Introduction 



The specific features of lepton mixing [TJ [2] inspire the search for lepton mass matrices 
which might reflect some of those features in a natural way. In our previous paper [3] 
we have constructed a model predicting trimaximal lepton mixing. This is defined by 
\U e 2\ 2 — [U^l 2 — \U T 2\ 2 = 1/3; where U is the lepton mixing matrix. Trimaximal mixing 
is a less stringent requirement than tri-bimaximal mixing [JJ, which furthermore imposes 
\Ufjzl = \U T z\ and U e3 = 0. In this paper we shall further constrain our model of [3] by 
adding to it the non-standard CP symmetry first introduced in [5] J]] 

Let us list the steps in the construction of the model and, simultaneously, illustrate 
how we proceed to reduce the number of parameters in the lepton sector: 

• Neutrinos are assumed to be of the Majorana type. The smallness of the neutrino 
masses is explained through the (type I) seesaw mechanism [7j. 

• Family symmetries — the discrete group A (27) together with three 7L>2 groups — are 
imposed, and afterwards broken softly and spontaneously in a sophisticated way. 
These symmetries justify specific forms for the charged-lepton mass matrix Mi, 
which is diagonal, and for the neutrino Dirac mass matrix Md, which is proportional 
to the unit matrix. The family symmetries also reduce the number of Yukawa 
couplings to a minimum and strongly constrain the Majorana mass matrix of the 
right-handed neutrinos Mr. Since Mi and Mp are both diagonal, lepton mixing 
originates solely in Mr. 

• We get rid of two phases by assuming invariance of the Lagrangian under the non- 
standard CP transformation of [5J, which is eventually broken spontaneously at the 
electroweak scale. 

The specific construction of the model will be explained in the next section; we anticipate 
that it produces the three-parameter light-neutrino mass matrix [3] 
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in the basis where Mi is diagonal. Our model thus predicts the neutrino masses and 
lepton mixing — a total of nine observables — in terms of just three parameters — the real 
numbers x, y, and z. The contribution of the non-standard CP transformation to the 
form of M v is to constrain these numbers to be real; in our original model [3] they could 
be complex and, therefore, the number of physical parameters in Ai v was five. 
With the M. v of equation ([!]), trimaximal mixing is realized since 
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1 Earlier, that symmetry had already been proposed in [B], but only for neutrinos and not within a full 
model. 
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where A = x+2z. Furthermore, as a consequence of the non-standard CP transformation, 
My fulfills [5] 

/ 1 \ 

SM U S = M* U , where 5= 1. (3) 

\0 1 0/ 

Let us discuss the predictions of equation (CQ) for the neutrino masses and for lepton 
mixing. We follow the convention in the Review of Particle Physics [H] for the param- 
eterization of the lepton mixing matrix. As shown in [5], it follows from equation ([3]) 
that 
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i.e. maximal atmospheric-neutrino mixing and maximal CP violation. Moreover, there 
are no CP- violating Majorana phases, i.e. the values of those phases are trivial (0 or it). 
On the other hand, it follows [3] from equation (T5]) that 
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which relates S\2 with s 13 and provides the lower limit 
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on Si2. The mass matrix does not determine the type of neutrino mass spectrum; 
it could be normal — wherein the smallest neutrino mass m s = mi — or inverted, i.e. 
m s — m 3- Still, equation ([T|) does fix the absolute neutrino mass scale as a function of 
both Aml tm 
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m{\ and sf 3 ; indeed, 
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(This is a result from [3] specialized to = 1/2.) In summary, due to M. v containing only 
three parameters, there are in our model six predictions for the nine physical observables 
following from M. v : 

• the two Majorana phases, the Dirac phase S, and the atmospheric mixing angle are 
all fixed; 

• the solar- neutrino mixing angle is a function of s 2 3 through equation ([6]); 

• the smallest mass m s is a function of both sf 3 and Am 2 tm through equation (jSJ). 

In that sense, the parameters x, y, and z in M. v can be traded for sf 3 , Am 2 tm , and 
Am 2 & = m\ — m\. 
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We want to stress that the M. v of equation ([T]) is different from the mass matrices 
based on \x-r interchange symmetry. With that symmetry one obtains S13 = 0, therefore 
CP violation in neutrino oscillations is absent]! In the model of this paper, on the other 
hand, S13 does not need to vanish and CP violation in neutrino oscillations is maximal. 
Recent fits to the oscillation data [2] indicate the possibility that sf 3 > at the low scale; 
this would, if confirmed, disfavour a /i-r-symmetric M. v . 

We also want to stress that the most general M. v satisfying equation (|2J) is 



( r + s u t ' 

u r + t s 
V t s r + u J 



with r,s,t,u& C 



(9) 



and A = r + s + t + u. The M. v of equation ([9]) contains seven parameters: the moduli of 
r, s, t, u and their relative phases. Our light-neutrino mass matrix is clearly a much 
restricted version of (jHJ). Notice that, with the M. v of equation (jHJ), tri-bimaximal lepton 
mixing is attained when t = u. 

This paper is organized as follows. In section [2] we discuss the model of [3] under the 
constraint of the non-standard CP symmetry of [5J, demonstrating that the neutrino mass 
matrix ([T]) is obtained. In section [3] we study the effect of the one-loop renormalizat ion- 
group evolution on Ai u and show that, in our model, the form ([TJ) of A4 U is not changed 
by that evolution, so that the ensuing predictions hold irrespective of the energy scale. 
The computation of the effective mass for neutrinoless f3(3 decay is the subject matter of 
section @J We present our conclusions in section [51 



2 The model 

The model that we have put forward in [3J is an extension of the Standard Model, with 
gauge group SU{2) x £7(1). The lepton multiplets are the standard left-handed SU{2) 
doublets D aL = (v a L, Oi) T , the right-handed charged-lepton 577(2) singlets c*r, and four 
right-handed-neutrino SU{2) singlets v a R, u 0R (a = e, //,r). The scalar sector of the 
Standard Model is also extended to four SU(2) doublets <p a , <fio, together with a complex 
gauge singlet S. 

Next we discuss the symmetries of the model. The two 3x3 matrices 
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satisfy F 3 = T 3 = 1 and do not commute. Together they generate one of the two three- 
dimensional faithful irreducible representations of the discrete group A (27) [10]; the other 
one is generated by F and T*. We display in tabled] the way in which the multiplets with 
index a transform under this horizontal symmetry A (27). The multiplets without index 

2 For a model with a three-parameter neutrino mass matrix based on /i-r interchange symmetry, see [9] . 
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Table 1: Transformation properties of the A (27) triplets under F and T. 

a are T-invariant and transform as 

vqr -> uv 0R , 0o ^0o, S -> cuS 1 under F. (11) 

Next we introduce three symmetries z ejAtjT as [H] 

z Q : ckh — > -a R , 4> a -> -0 a , (12) 

and all other multiplets remain unchanged. Finally, we come to the non-standard CP 
transformation [5], defined by 

D aL -> iS a p-f CDp L , a R -> iS a pj C[3 R , v aR -> iS a pj Q Ci?J R , u 0R -> ij Ci?Q R , (13) 

for the fermions (C is the charge-conjugation matrix) and 

0q — > 5 Q/3 0^, 00 ^0o, 5 — > S* (14) 

for the scalars. The matrix 5 is given in equation (13"1) . 

The multiplets and symmetries lead to the Yukawa couplings [3] 

£y = -yi ^ D aL a R (j) a - y A DaL^aR (*t 2 0o) + y ^orC -1 ^/? 5 + H.c. (15) 

In the Majorana mass terms we allow soft breaking of the symmetry T, but not of F. 
Therefore, [3] 

Note the consequences of the non-standard CP symmetry: 

J/i,3/4,!te,M 0) M ll M 2 eR. (17) 

Since the first term in £y has one common coupling constant yi, upon spontaneous 
symmetry breaking one needs three different vacuum expectation values (VEVs) v a = 
(0°) ^° ^ e a ^^ e ^° accoun t f° r the different charged- lepton masses: 

m e : m^: m T = \v e \ : : |u r | . (18) 
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From equations (1151) and (I16p . we find the Majorana and Dirac neutrino mass matrices 
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(19) 



respectively. We have defined = y$v* s and a = ?/4t>o, where vs is the VEV of S and t>o 
is the VEV of the neutral component of 0o- 

With M v = -M^M^Mn, the mass matrix ([!]) is obtained. Since M 0i i i2 e E, 
it follows that the parameters x and z of M. v are real as well — see [3]. On the other 
hand, ?/ oc 1/Mjv [3]; therefore, the question whether y is real boils down to the question 
whether the minimum of the scalar potential V occurs for a real v$- For this reason, we 
next investigate V. 

The terms of dimension four in V must be invariant under F, T, and z EjAt)T . Terms of 
dimension three are invariant under F and z e At T , but they may violate T, cf. equation (TT6l) . 
If one requires terms of dimension two to be invariant under z e>fl)T as well, then there are 
four U(l) symmetries in the potential, one for each Higgs doublet; these symmetries are 
all spontaneously broken, resulting in three physical Goldstone bosons. To avoid this 
problem we break the z a softly by terms of dimension two. As for the symmetry F, we 
want it to be spontaneously broken through different VEVs v e ^ T \ in order to guarantee 
that this happens, we must admit F to be softly broken by terms of dimension two too. 
The only symmetry respected by the full scalar potential is the CP symmetry. Thus, 
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A 7 {4>\<p e + <f>\ 4>„ + 40 r ) + A 8 0j0o 
+A 9 [S 2 (0t0 e + u 2 ^, + cu0t0 T ) + S 2 * (0l0 e + cu^ + ^ 2 0t0r)] 

+fl [S (0l0 e + + ^ 2 0t0r) + S* (0l0 e + UJ 2 ^ + UJ4>\<P T )] 

+H \S\ 2 + A \S\ 4 + fn (S 3 + S* 3 ) + /i 2 (S 2 + S* 2 ) + /i 3 (S + S*) 

i,j=0,e,fj,,T 



(20) 



All the parameters in V are real; the only exception are the parameters of the 4x4 matrix 
M.Ay, which has the structure 



(a b c c* \ 

b d e e* 

c* e* f g 

\ c e g* f ) 



with a, b,d, f G R, 



(21) 
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in the basis (0o, 4> e , <P^, <j> T )- Note that c, e, and g are in general complex. These conditions 
derive from hermiticity together with CP invariance. Note that M.^ is still general enough 
to allow for different VEVs Vi (i = 0, e, /i, r). 

Since the VEVs Vi are of the electroweak scale, whereas vs is of the seesaw scale, a 
fine-tuning is necessary in V, with extremely small coupling constants A7, As, A9, and \x. 
This is an unpleasant feature of the model, which is shared, though, by all other models 
with two widely different mass scales. 

Now we address the conditions for obtaining a real VEV of S. For this purpose we 
only have to consider the penultimate line of equation ( J20l) . Obviously, if 



sign Hi = sign /i 3 and /x 2 < 0, 



'.22) 



then the absolute minimum of V is attained for real v$ (vg > for fii < and v$ < for 
Mi > 0). 

We find it useful to summarize the symmetry breaking in our model: 

• At the seesaw scale, T is softly broken by terms of dimension three in Lyia], while 
F is spontaneously broken by the VEV of 

• At the electroweak scale, T, F, and the z a are all softly broken in the matrix M.^. All 
symmetries, including the non-standard CP symmetry, are broken spontaneously 
at this scale. 

The CP symmetry is broken by 7^ |u T |; the different VEVs are necessary in our model 
for having m M 7^ m T — see equation ( 1T81) . Now we want to show that indeed there is no 
violation of CP for degenerate masses m M = m T . In that case we may transform 



M R -> KMrK with K = 

where K' is a 2 x 2 unitary matrix. Choosing 
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one sees that 
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(25) 



becomes real. Thus, it is indeed \v^\ 7^ \v T \ which breaks the CP symmetry, at the 
electroweak scale. Notice that the phases of the VEVs are irrelevant in the breaking of 
CP. 



3 Notice that the VEV of S is crucial in our model in order to obtain Mjy 7^ 0. If it were not for 
Mn 7^ 0, the seesaw mechanism would be unable to suppress the masses of all three light neutrinos, 
because det Mr oc Mn- 
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3 The renormalization-group evolution 



Throughout this section, we use indices = 0,e,fi,r to refer to the four Higgs 

doublets of our model. 

We use the formalism of our paper [12]. In that paper we have considered a multi- 
Higgs-doublet extension of the Standard Model supplemented by effective dimension-five 
operators 

o*= E - °£*t) c- 1 ow; - m$) , (26) 

a,f3=e,/j,,T 

where the k^) are coefficients with dimension —1. This is sufficient for our purposes 
despite the occurrence, in our model, of a scalar singlet S. The reason is that S is 
integrated out at the seesaw scale iur, because S has VEV and mass of order tur. 

In our model there are symmetries F, T, and z J As emphasized at the end of the last 
section, both F and T are broken at the seesaw scale thr, while the three symmetries z a 
stay valid below rriR and are broken, both softly and spontaneously, only at the electroweak 
scale mp. Thus, at the relevant energy scales, i.e. in between itlr and mp, the symmetries 
z a hold. We therefrom conclude that in our model the only non-zero matrices k^' are 
those with i — j. 

Upon spontaneous symmetry breaking the light-neutrino mass terms are 

~ E ( M »U "Lc-Wp with ^ = ;£ w «a (27) 

OL,P i,j 

where V{ is the VEV of 0°. Since only the are non-zero, we have, at any energy scale, 

M„ = 2j2 v iK (il) - (28) 

i 

The matrices are symmetric, as is obvious from equation ff26l) . 

At the high scale m R , only the Higgs doublet 4> (and also the scalar singlet S) has 
Yukawa couplings to the right-handed neutrinos, while the other three Higgs doublets 4> a 
have no such couplings. We therefrom conclude that, at the scale tjir, 



K W = ^ ^) = 03 X3 fori = e ,/,,T. (29) 
2^o 



This is the initial condition for the renormalization-group (RG) running, with M. v given 
by equation flTJ). 

Let us write the scalar potential of the Higgs doublets in the form 

V = quadratic terms + \jki4 > \4>j 4>\4>i- (30) 

i,j,k,l 



4 There is also the CP symmetry. The reasoning and the conclusions in this section are, however, in- 
dependent of the presence or absence of that CP symmetry. Thus, the present section and its conclusions 



apply as well to the model of [3] . 
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Let us furthermore write the Yukawa couplings of the Higgs doublets to the right-handed 
charged leptons in the form 

£ ^ = - 2 2 AxL/Mi + H.c. (31) 



Then, the RG equations for the evolution of the matrices kS u ' are 

i67r 2^_ = (_3^ + 2Ti .^(u) + fi .(n)p + pT ft .(n)_ 2 pu), 

(32) 



Here, t = In /x is the logarithm of the mass scale, g is the SU(2) gauge coupling constant, 
and 



T ifc = tr^ T y fc J, (33) 
P = \Y. Y i Y l ( 34 ) 

3 

If we define the 3x3 matrices 

Po = 3X 3, P e = diag (1,0,0), P M = diag (0,1,0), P T = diag (0, 0, 1) , (35) 
then we see in equation (115]) that, in our model, at the scale uir, 

Yi = 2/i Pj for z = 0, e, /i, r, (36) 

and, therefore, 

I ,2 

P ee = P w = Prr = M 2 and all other Ty = 0, P = -i— l 3x3 . (37) 

Now one can easily check that the matrices Yi remain of the form (1361) at all energies 
when they evolve with the RG equations for the Yukawa couplings. We therefore have 

= (-3/+| 2/l | 2 )^°) + 4^A^ fe ), (38) 

3 

A K {aa) 

16vr 2 ^^ = (-3g 2 + 3\y 1 \ 2 )K^-2\y 1 \^K^\P a }+Aj2^ ja ^\ (39) 

3 

where the anti-commutator of matrices A and B is denoted {A, B}. 

We observe that the coefficients \ i 3 i are particularly important in equations (1381) 
and (139]) . According to [12], the RG equation for those coefficients is 

167T 2 = 4 £ (2A iiW A^ + AywAaw + Kkij^mj + W*y + A^y) 

- (% 2 + 3</' 2 - 2T« - 22^) Aijij (40) 



whenever % ^ j. In equation (I40p . g' is the U(l) gauge coupling constant. A little 
contemplation of equation (HOI) allows one to ascertain that, if all the A^ with i ^ j 
vanish at the scale m^, then no such non- vanishing coefficient will ever arise through the 
RG evolution in our model. 

Now, this is precisely the situation that occurs. As seen in equation ( 1201) . at the scale 

i.e. at the initial condition for the RG evolution, there is no term (^4>\(pj^ with i ^ j 

in the scalar potential. Therefore, all coefficients A^ with i ^ j vanish at the scale 
/i = m R , and, through equation (HOI) , at all other scales too. 

Equations (1381) and (|39|) . together with the initial condition (1291) . then allow us to 
ascertain that, in our model, all matrices k^> except vanish at all energy scales. The 
sole non-vanishing matrix evolves according to 

167T 2 = (-3/ + M 2 + 4A ooo) k (00) , (41) 

and therefore it does not change its form along the RG evolution. This means that all 
the predictions of our model are RG-invariant. 

In summary this result comes from the fact that the symmetries of our model suitably 
constrain both the Yukawa and the quartic Higgs couplings, and that soft and spontaneous 
breaking at the seesaw scale — which is necessary in our model for obtaining the desired 
M. v — has no impact on their RG equations. 



4 The effective mass in neutrinoless (3(3 decay 

In our model the effective mass relevant for neutrinoless (3(3 decay is 

mpp = \(M u ) ee \ = \x + y\. 
In order to compute this we have to diagonalize M. v . We define [1] 



^hps = 



2/Vq i/Vz o \ 

-1/V6 1/V3 -l/y/2 

l/y/3 1/V2 J 



and compute 



/ p + q 



Therefore, defining 



iq 



\ 



x + 2z 
\ iq p — q J 



with p = x — z and q 



3y 



V = Ukps Di K, with A = diag (1, 1, i) and K 



( cos a sin a ^ 

1 
y — sin a cos a J 



(42) 



(43) 



(44) 



(45) 
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it is evident that the mass matrix M. v can be diagonalized as 

V T M u V = dmg(X 1 ,X 2 ,X 3 ), (46) 

where the neutrino masses are given by rrij = \Xj\ (J = 1, 2, 3) and the full diagonalization 
matrix is 

U = Vdiag (7^,7/2, 773) , (47) 

where r)j = 1 for Xj > and rjj = i for Xj < 0. 
Then, 

(M0« = Al (K*l) 2 + A 2 (K*2) 2 + A 3 (Ks) 2 

2Ai + A 2 , . . N 2 sin 2 a 

- —3— + (A3 - AJ —3— 

2Ai + A 2 , 1 - cos 2a 

- — g— + (A3-A0 

Ai + A 2 + A 3 (Ai - A 3 ) cos 2a 

3 3 1 J 

Since 

/ cos a — sin a \ / p + a —q \ ( cos a sin a , 

diag(A 1 ,A 3 )= . 1 - > ( 49 ) 

sin a cos a / \ — q —p + q I \ — sina cos a 



one finds that 



sin 2a = — 9 (50) 



ca/P 2 + g 2 

cos 2a = —===, (51) 



where e = ±1. Thus, 
hence 



Ai = q - ey/p* + g 2 , (52) 
A3 = g + eyV + g 2 , (53) 

A1A3 = V < 0, (54) 



V = W-A1A3, (55) 
where r\ = ±1. From equations (|5Tj) (T53T) . we obtain 

cos2a = - ^— . (56) 

A3 — Ai 

Returning to equation (jl8|) . one obtains 
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Now we must examine the matter of sign ambiguities in (Ai u ) ee . The three Xj may 
be either positive or negative, but they are subject to the condition A1A3 < 0. Thus, 

Ai = (mi, A 2 = em 2 , A 3 = — (m 3 , with £ = ±1 and e = ±1. (58) 

On the other hand, 

Ity 2 = (59) 
is experimentally known to be very small, or even zero. Hence, 

cos 2a = 1 - 3 \U e3 \ 2 > 0. (60) 

Therefore, from equations (J55|) and (J55|) . 



or 



Therefore, 



~ 2 :'™ > o, (6i) 

A3 — Ai 

2?]v /= AiA3 



C (mi + m 3 ) 



> 0. (62) 



^a/-AiA 3 = CV m i m 3- (63) 

Thus, from equation (l5"Tj) . 

3 (^f)ee = £m 2 + C (mi - m 3 + 2 v /mTm 3 ") . (64) 
There are therefore two possibilities: 



either mpp = - \m 2 + mi — m 3 + 2y/mims\ , 
3 

1 . , . 

or mpp = - |m 2 — mi + m 3 — 2^/m\m 3 \ ■ 



(65) 



Both these possibilities exist independently of whether the neutrino mass spectrum is 
normal or inverted. Thus, besides fl65l) . 



atm; 



either m\ = m s , m 2 = \/m 2 s + Am|,, m 3 = a/ m 2 + Ar 

(66) 

or mi = ^/m 2 + Am 2 tm , m 2 = a/ib 2 + Am 2 tm + Am|, m 3 = m fl . 

In figure [1] we have plotted m/3/3 as a function of s 2 3 , taking into account all four 
possibilities explained above. The input values for the mass-squared differences are 
Am 2 tm = 2.40 x 10" 3 eV 2 and Awig = 7.65 x 10 -5 eV 2 , the mean values given in the 
second paper of [2]. 
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Figure 1: The effective mass in neutrinoless (5(5 decay, m^, as a function of \U e z\ ■ The 
full lines and the dashed-dotted lines correspond to a normal and an inverted neutrino 
mass spectrum, respectively. We have fixed the neutrino mass-squared differences at their 
mean values given in the second paper of [2J. 

5 Conclusions 

We have presented in this paper a model whose symmetries force the charged-lepton mass 
matrix to be diagonal while generating the highly predictive three-parameter neutrino 
mass matrix of ([I]). It follows from this neutrino mass matrix that the lepton mix- 
ing matrix U has only one free parameter, which may be chosen to be S13, while both 
atmospheric-neutrino mixing and CP violation are fixed and maximal. Note that, while 
U displays a maximal Dirac phase, it has vanishing Majorana phases. Since the model 
leads to trimaximal mixing, s\ 2 must be larger than 1/3 — see equations (jSJ) and (|7j). This 
is slightly disfavoured by the present data, but the value 1/3 is still within the 2a range 
for s\ 2 - In any case, the correlation between s\ 2 and is a crucial test of trimaximal 
mixing. 

In our model, the symmetry which leads to maximal atmospheric-neutrino mixing and 
maximal CP violation is the non-standard CP transformation given by equations ( fTBl 
and dHJ). A CP symmetry of this type has the curious property that m M 7^ m T is an 
effect of its spontaneous breaking, as was noticed earlier in [5], IT3] . 

Given s\ 3 and Am^ tm , the mass matrix (CQ) fixes the absolute neutrino mass scale — see 
equation (jBJ) — but does not determine the type of neutrino mass spectrum. Because of 
a sign ambiguity there are, for each type of spectrum, two possibilities for the effective 
mass mpp of neutrinoless (5(5 decay as a function of s\ z . If the neutrino mass spectrum 
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is of the inverted type, one of the possibilities for mpp is within the projected range of 
future experiments. 

We have also presented in this paper a mechanism for the one-loop renormalizat ion- 
group stability of the neutrino mass matrix. Indeed, although some soft and spontaneous 
breaking of symmetries occurs already at the seesaw scale, the symmetries of our model 
are such that, in between the seesaw and the electroweak energy scales, only the effective 
dimension-5 neutrino-mass operator associated with the Higgs doublet 0o is non-zero, 
but that Higgs doublet is different from the Higgs doublets 4> e ,^,T which give mass to the 
charged leptons. Therefore, the form of the neutrino mass matrix is not RG-distorted by 
the fact that all three charged leptons have different masses. 
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